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I. INTRODUCTION 



Fluid dynamical models are widely used to describe ultra-relativistic heavy ion collisions. Their advan- 
tage is that one can vary flexibly the Equation of State (EoS) of the matter and test its consequences on 
the reaction dynamics and the outcome. For example, the only models which may handle the supercooled 
QGP are hydrodynamical models with corresponding EoS. In energetic collisions of large heavy ions, es- 
pecially if a Quark-Gluon Plasma (QGP) is formed in the collision, one-fluid dynamics is a valid and 
good description for the intermediate stages of the reaction. Here, interactions are strong and frequent, 
so that other models, (e.g. transport models, string models, etc., that assume binary collisions, with 
free propagation of constituents between collisions) have limited validity. On the other hand, the initial 
and final, Freeze-Out (FO), stages of the reaction are outside the domain of applicability of the fluid 
dynamical model. 

Thus, the realistic, and detailed description of an energetic heavy ion reaction requires a Multi Module 
Model, where the different stages of the reaction are each described with a suitable theoretical approach. 
It is important that these Modules are coupled to each other correctly: on the interface, which is a three 
dimensional hyper-surface in space-time with normal dcr^, all conservation laws should be satisfied, e.g. 
[T^duu] = (here the square brackets means difference between new and old phases or modules), and 
entropy should not decrease, [S^da^] > 0. These matching conditions were worked out and studied for 
the matching at FO in detail in Refs. (1]-^. 

We would like to discuss the entropy condition in more detail. Obviously, the number of degrees of 
freedom and correspondingly the entropy density is reduced during hadronization process. So, how can 
we avoid decreasing of the entropy? Two scenarios have been proposed. The first one is the gradual 
hadronization scenario i.e., the hadronization is so slow that during this process the volume of the system 
becomes considerably larger to compensate for the reduction of entropy density. If this would be so, 
our long living, gradually expanding QGP should be observed in HBT experiments, e.g., as a peak in 
the Rout/Rside ratio jj]. The preliminary data from STAR and PHENIX do not support this scenario 
|J. The second possibility is the fast hadronization from supercooled QGP This hypothesis can be 
checked only in hydrodynamical models, which use the EoS as direct input. 



After hadronization and FO, matter is already dilute and can be described well with kinetic models. 

The initial stages are more problematic. Frequently, two or three fluid models are used to remedy the 
difficulties and to model the process of QGP formation and thermalization |]Io|-|l^|. Here the problem is 
transferred to the determination of drag-, friction- and transfer- terms among the fluid components, and 
a new problem is introduced with the (unjustified) use of an EoS in each component in a nonequilibrated 
situations where an EoS is not defined. Strictly speaking this approach can only be justified for mixtures 
of noninteracting ideal gas components. Similarly, the use of transport theoretical approaches assuming 
dilute gases with binary interactions is questionable, as, due to the extreme Lorentz contraction in the 
Center of Mass (CM) frame, enormous particle and energy densities with the immediate formation of a 
perturbative vacuum should be handled. Even in most parton cascade models these initial stages of the 
dynamics are just assumed in form of some initial condition, with little justification. 

Our goal in the present work is to construct a model, based on the recent experiences gained in string 
Monte Carlo models and in parton cascades. One important conclusion of heavy ion research in the last 
decade is that standard 'hadronic' string models fail to describe heavy ion experiments. 

All string models had to introduce new, energetic objects: string ropes |]l3| , |l4| , quark clusters [fj"5| , 
or fused strings [ fl6| , in order to describe the abundant formation of massive particles like strange an- 
tibaryons. Based on this, we describe the initial moments of the reaction in the framework of classical 
(or coherent) Yang-Mills theory, following Ref. [^0) assuming a larger field strength (string tension) than 
in ordinary hadron-hadron collisions. For example, calculations both in the Quark Gluon String Model 
(QGSM) Jl7|-p9[ and in the Monte Carlo string fusion model [[l6| indicate that the energy density of 
strings reaches 8 — 10 GeV/ fm already in SPS reactions, nearly 10 times more than the tension used in 
standard, 'hadronic', string models where a ~ 1 GeV/ fm. In addition we now satisfy all conservation laws 
exactly, while in Ref. [ppf infinite projectile energy was assumed, and so, overall energy and momentum 
conservation was irrelevant. Thus, in this approach for the first time the initial transparency/stopping 
and energy deposited into strings and "string ropes" will be determined consistently with each other. 
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Recent parton kinetic models pl|j2^ ] indicate that quark and gluon density saturation takes place in a 
very short time - r sa t = 0.09 — 0.27 fm/c for LHC - SPS energies [P2| , and pressure needs t p = 1 — 5 fm/c 
to be established . More importantly the first experiments at RHIC yield strong elliptic flow, which 
cannot be reproduced in any other model, except in fluid dynamical models with QGP EoS |^3| . This is a 
strong experimental indication that transverse pressure builds up early in these reactions, in a few fm/c, 
and strong stopping is also necessary to create strong flow before Freeze Out, which usually happens 
when the system size is not more than 10 fm. We present initial conditions for i/ a {, = 2 — 5 fm/c, which 
is in agreement with previous estimations as well as with data. 

We do not solve simultaneously the kinetic problem leading to parton equilibration, but assume that 
the arising friction is such that the heavy ion system will be an overdamped oscillator, i.e., yo-yoing of 
the two heavy ions will not occur, as all recent string and parton cascade results indicate. 

II. FORMULATION OF MODEL 

Our basic idea is to generalize the model developed in [^o| for collisions of two heavy ions and improve 
it by strictly satisfying conservation laws p^-|26|]. First of all, we would create a grid in the [x,y] plane 
(z is the beam axes, [z,x] is the reaction plane). We will describe the nucleus- nucleus collision in terms 
of steak- by-streak collisions, corresponding to the same transverse coordinates, {xi,yj}. We assume that 
baryon recoil for both target and projectile arise from the acceleration of partons in an effective field 
F^ v , produced in the interaction. Of course, the physical picture behind this model should be based on 
chromoelectric flux tube or string models, but for our purpose we consider as an effective Abelian 
field. The most important consequence of the non- Abelian fields, i.e., its self interaction and the resulting 
flux tubes of constant cross section, are, nevertheless, reflected in our model: assuming that the field is 
one-dimensional. The fields generated by the colliding streaks are of constant cross section during the 
whole evolution, and only their lengths increase with time. As the string tension is constant, the energy 
of the string increases linearly with its increasing length. The single phcnomcnological parameter we use 
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to describe this field must be fixed from comparison with experimental data. 
We describe the streak-streak collision using conservation laws: 



i i 

^E^ = - z = 1 ' 2 ' ( 2 ) 

i 

where N]* is the baryon current of ith nucleus, qt is the color charge and it will be discussed in more 
detail later. We are working in the Center of Rapidity Frame (CRF) , which is the same for all streaks. 
The concept of using target and projectile reference frames has no advantage any more. We will use the 
parameterization: 

ATf = riiuf , = (coshj/i, sinhy;) . (3) 

T^ u is the energy-momentum flux tensor. It consists of five parts, corresponding to both nuclei and free 
field energy (also divided into two parts) and one term defining the QGP perturbative vacuum. 

T^ = E^+^t = Eh(( 1 + c o)<<- c o^1+ T F,i] +B(T , i = l,2. (4) 

i i 

Here B is the bag constant, the equation of state is Pi — c^ei, where and Pi are energy density and 
pressure of QGP. 

Within each streak we form only one flux tube with a uniform field strength or field tension, a, from 
the target to the projectile. For practical purposes we, however, divide this field into two spatial domains, 
a target and a projectile domain, (i = 1,2), separated at a fixed point, z sep , so that a\ = <r 2 = a. The 
choice if this point will be specified later. (The field is constant and the only change is that it extends 
with time at its two ends.) 

In complete analogy to electro-magnetic field 

Fr = d»A»-d»A?=(° , (5) 
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*i = d 3 A° - d°A* , (6) 

£f,< = -jFi^vF^ v . (8) 

In our case the string tensions, Uj, will be constant in the space-time region after string creation and 
before string decay. The creation of fields will be discussed later in more detail. 
To get the analytic solutions of the above equations, we use lightcone variables 

(z,t) -» (x + ,x~), x ± =t±z. (9) 

Following pp[|, we insist that ei,$/i,7ii are functions of x~ only and e2,j/2,"2 depend on x + only. 
In terms of lightcone variables: 

N± = N hT = ni(u? ± uf) = me*** , (10) 

t - ( T ? + T ?~ \ - 1 ( hi+e2vi hi ~ Ut mi 

where 

h i+ = (1 + c^e* , = (1 - cg)e, . (12) 

The other tensors in the light cone variables are: 

t - = (bo)' (14 » 

The energy-momentum tensor for free field in the light cone variables is: 



At the time of the first touch of two streaks, t = 0, there is no string tension. We assume that strings are 
created, i.e., the sting tension achieves the value a at time t = to, corresponding to complete penetration 
of streaks through each other (see Fig. |]). 

Streak-streak collision 




FIG. 1. Streak-streak collision, t = at the time of first touch of streaks, t = to corresponds to complete 
penetration of streaks through each other. At this time strings are completely created, i.e., string tension reaches 
an absolute value a = A (f^-) 2 noVhh (|26|). 



III. CONSERVATION LAWS AND STRING CREATION 



In lightcone variables eq. (g) may be rewritten as 



cLiVF + d+N+ = 



(16) 



So, we have a sum of two terms, each depending on different independent variables, and the solution can 
be found in the following way: 



d+N+ = -a , 



cLiVr = a, 
N- = ax- + (JVf ) , ^ = -ax+ + (N+) 



(17) 



G 



where the index )o indicates the initial proper density, which is the normal nuclear density: no = 
0.145 /m~ 3 . Since both iV-f and are positive (and also more or less symmetric) we can conclude 
that for our case a — 0. 
Finally 

Ni = me- yi = n e yo , 7V 2 + = n 2 e V2 = n e Vo , (18) 

rii = n e Vo+Vl , n 2 = n e v °- y2 . (19) 

where yo (—yo) is the initial rapidity of nucleus 2 (1) in the center of rapidity frame (CRF), respectively. 
The other components are given by eq. (|l0|). 

Let us make analogy to electro-magnetic field, where two charges q\ and —q 2 move in the opposite 
directions, creating string- like field between them, E — (0, 0,1?), which is constrained transversally into 
a constant cross section. The Z-axis goes through charges q\ and —q 2 and directed from q\ to —q 2 (let us 
assume that we have such a field configuration). So, forces acting on our charges, q\E and —q 2 E, have 
opposite sign and both are working against the expansion of the " string" . In our effective model wc use 
color charges, and assume that the vectors of these color charges point in the opposite directions in the 
color space pq , so that the forces acting on both target and projectile partons are opposite both stopping 
the expansion of the streak. As our field strength (string tension, a) is not yet defined we normalize the 
charges to unity: 

qi = —q 2 = 1 > , while <j\ = <r 2 = a . (20) 

Then we have the forces acting in z direction: q\U\ — a , and q 2 a 2 = —a. Notice again that after string 
creation fields <Ji(x) and a 2 (x) are spatially separated as well as the baryon densities, n\ and n 2 \ i.e., 
after complete penetration of the initial streaks through each other (see Fig. 1), a 2 acts on the partons 
on the right side of the separating point z sep = (h — l 2 )/2, while the <j\ acts on those on the left side. 
In the absence of matter, in the middle both fields are identical, so the exact position of the separating 



point does not play any role until it does not enter the target or projectile matter. The fields o\ and 
o 2 are generated by the corresponding 4-potentials Ai , which are different and spatially separated in the 
same way. 

As it was described above we do not generate the chromo-electric field self consistently, as a product 
of color currents, which a affected also by the field. Our effective fields are external with respect to 
colliding partons, that why we can use the expression ([f5|) for the field energy. On the other hand, if we 
want to satisfy the conservation laws, we must generate our effective fields in the collision transferring 
energy from matter to field. It's possible to define new conserved quantities based on eq. Using the 
definition of F^ v , eq. (||), we can rewrite eq. (jlj) as 

d ^ v = J2 FfqiNi^ = £ © ^ ( A i N ^) ~ - d v «iV t> ) + AS&N^] . (21) 



The solutions for N 1 and , eq. (18), show that the second term vanishes. The fourth term is a 



vector —{A l d-N^ , j4;j~(9 + iV 2 ) in lightcone coordinates. So, if we impose the conditions 

A^=0, A+=0 (22) 
we can define a new energy-momentum tensor T Miy , such that 

fyT"" - , (23) 

fnv = ^ fr + T p Z t = Pi" - 1i A i N i + g^liAfNia) + Bg»" (24) 

i i 

To satisfy the above choice of the fields, (|2C|), and imposed conditions ( p2[ ) together with the Lorentz 
gauge condition, d°A® — d 3 Af = or d~A^ — d + Af = 0, we take the vector potentials in the following 
form: 

Ai = 0, A^ = = -<rx+, ^ 

A% = <J2X~ = ox~ , A\ = . 



Notice that the above choice differs from the one which was initially proposed in Refs. [[24|-[26| , what 
causes the changes in the expressions related to field creation in this section - eqs. (|27],|3l],|3^) - but will 
not affect the analytic solution of the model - eqs. (134 
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In our calculations we used the parameterization: 

<T = A(^) 2 n ^hh, (26) 

where m is the nucleon mass, and l\ and 1% are the initial streak lengths (see Fig. |l|). We are working 
in the system, where Ti = c = 1, so a has a dimension of length^ 2 = energy /length. The typical values 
of dimcnsionlcss parameter A are around 0.06 — 0.08. Notice, that there is only one free parameter in 
parameterization (p6|). The typical values of a are 4 — 10 GeV/fm for Eq = 65 GeV per nucleon, and 
a pa 6 — 15 GeV I fm for £o = 100 GeV per nucleon. These values are consistent with the energy density 
in non-hadronized strings, or "latent energy density" which is on the average 9 GeV/ fm 3 [pT|-^9[. 

Using the exact definition of Af, eqs. (p5l), eqs. (pd|,^4],^5|,p0|j2^) and transformation matrixes from 



Appendix C we obtain 



f++ f+- \ _ 1 / h l+ e 2 y^ tn- 

f-+ f— j ~ 2 1 hi- h 1+ e~ 2 y i 



1 / h 2 +e 2y2 h 2 - \ If a 2 2B 

2 I ha- h 2+ e- 2 ^ ) + 2 { 2B a 2 



(27) 



/ -(JX+N+ 0\ ( ax'N+ \ 
+ y crx+N^ J + I -ax-N^ J 

2 

Notice that perturbative term, £?, and free field energy, cover all the interacting volume, while energy 
densities of matter and baryon currents are separated in space. We also want to stress factor 1/2 in front 
of all terms in (it has been canceled by 2 near a in last two terms) - this factor was missed in Refs. 
|2~f| |2G| ] as well as in ^0|, although it does not affect the result since equations of motion, d^T^ v = 0, can 
be multiplied by any coefficient. The reason for it is a form of transformation matrixes between (t, z) 



and (+, — ) coordinates, which are presented in Appendix C 
Now the new conserved quantities are 



Qo = / f 00 dV = Ax A y / f 00 dz , (28) 
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!3 



T U W = Ax A y / T' 3 dz , (29) 



where the volume integral runs over the lengths of the both streaks and Ax A y is the cross section of 
the streaks. Notice that in the absence of the fields, before string creation and after string decay, the 
(Q°, Q 3 ) come back to (P°, P 3 ) - components of the 4-momenta of the system. 
We can rewrite energy-momentum tensor in (t, z) coordinates: 



foo f03 \ / ( ei -i-P^cosh 2 ^ - Pi (ei + Pi) cosh yi sinh yi 
T 30 f 33 / _ I (ei +Pi)coshyisinhyi (ei + Pi) sinh 2 y x + Pi 



(e 2 + P 2 ) cosh 2 y 2 - P 2 (e 2 + P 2 ) cosh y 2 sinh y 2 \ ( %■ + B \ 
(e 2 + P 2 ) cosh y 2 sinh j/i (e 2 + P 2 ) sinh 2 y 2 + P 2 J ^ ^ - B J 

ax+ ( N±-N+ Af-iV+ \ ax- ( N+ - N% - (N 2 ~ — N 2 ) 



(30) 



2 y-(N- +N+) -(N- +N+) J 2 yN+ + N- -(N++N-) J 

Based on conservation of Qo, Q3 we can calculate energy densities, ei(to), e 2 (io), at the moment t = to, 
when the string with tension a is created. These new quantities are used as initial conditions for our 
differential eqs. (|l|, ||). As shown in Appendix A - 



_ n m ^ 2 r + B l x +l 2 _ crn e vo 

ei(<o) = i +c r ( -) 2 (i^) - 4fe) 2 (1+c ^- < 31 > 

n m — +B h+li an e v " 

e2(o, = ^~er^^~*a)W 2 ' ( 1 

Here the ei(io) is a proper energy density at the time to, £0 is the initial energy per nucleon. We assumed 
transparency, i.e., that complete penetration happened so fast, that the fields, created during this time, 
did not have time to stop partons. So, the rapidities are 2/i( 2 )(io) = —yo(yo) correspondingly, and the 
proper baryon densities did not change. 

For x^ > xo we should solve eqs. (p3[), with boundary conditions 
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Ni(x~ =x ) = n a e Tyo N 2 ± (x+ = x ) = n a e 



±yo 



ht+(x = x ) =ei(t )(l + <^) h 2+ (x+ = x Q ) = e 2 (to)(l + c 2 ,) 



Vx{x =x Q ) = -y 



y 2 {x + =x ) = y 



(33) 



0\{x = xq) = a &2(x + = xq) = & 

qi{x~ = x ) = 1 q2(x + = x ) = -1 , 

where xo — 2to — \z(0)\ defines the string creation surface t — to, for parton or cell element in the position 
z = z(0) at the time t = 0. 

Let us present the complete analytical solution in the following form (for detailed calculations see 
Appendix B| ) 



(34) 



X — Xq 



(35) 



rii = uqc J e ' y 



(36) 



where x 1 = x , x 2 = x + , i, j — 1, 2 , i ^ j, and the notations are from Appendix B (B.7.B.9, B.13 - B.15 ). 
Then the trajectories of partons (or cell elements) for both nuclei are given by: 

x+{x-) = \z(0)\ + Q dx e 2 ^*) = 

(37) 



\z \-^( x -- Xo ) + ^+e- 2 ^ 



aa2 



2an e y 



1 



x 2 {x+) = |z(0)| + /* o e-^W = 
| Z( 0)| -£(*+- x ) + (£ + e"^) r 25 ^ (l - - ^a) 



(38) 



for parton or cell element in the position z = z(0) at the time t = 0. 
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IV. RECREATION OF MATTER 



If we let partons (or cell domains) evolve according to the above trajectories, eqs. (|37], [58|), they will 
keep going in the initial direction up to the time t = U t t urn , then they will turn and go backwards until 
the two streaks again penetrate through each other and a new oscillation will start. Such a motion is 
analogous to the "Yo-Yo" motion in the string models. Of course, it is difficult to believe that such a 
process would really happen in heavy ion collisions, because of string decays, string-string interactions, 
interaction between streaks and other reasons, which would be difficult to take into account. To be 
realistic we should stop the motion described by eqs. [58]) at some moment before the projectile and 
target cross again. 

We assume that the final result of collisions of two streaks, after stopping the string's expansion and 
after its decay, is one streak of the length AZj with homogeneous energy density distribution, e/, and 
baryon charge distribution, n/, moving like one object with rapidity yf. We assume that this is due 
to string-string interactions and string decays. As was mentioned above the typical values of the string 
tension, a, are of the order of 10 GeV/ fm, and these may be treated as several parallel strings. The 
string-string interaction will produce a kind of "string rope" between our two streaks, which is responsible 
for final energy density and baryon charge distributions. For simplicity we assume homogeneous baryon 
charge distribution. Notice, that in this way, after the decay of our "string rope" charges do not remain 
at the ends of the final streak, as it would be if we assume full transparency. The real situation may 
be more complicated: when the energy accumulated in the strong color fields will be finally released 
in a production of qq pairs and gluons, this process may noticeably change composition of matter as 
compared to the chemical equilibrium case |27j] . Therefore, matter created after the mutual stopping 
of interpenetrating streaks can not, in general, be described by the equilibrium EoS. The homogeneous 
distributions are the simplest assumptions, which may be modified later based on experimental data. Its 
advantage is a simple expression for e/, n/, yf. The first experimental results from RHIC do not show 
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transparency, rather most particle multiplicities as well as the elliptic flow show strong stopping and a 
peak around mid rapidity J23|. Furthermore, we describe the initial state, which is not directly observable 
in experiments, and a flat initial rapidity distribution may end up in both a forward-backward peaked 
and a centrally peaked distributions depending on several other circumstances. 

The final energy density, baryon density and rapidity, e/, n/ and yj, should be determined from 
conservation laws. The assumptions we made above oversimplify the situation, and do not allow us to 
satisfy exactly all conservation laws. The reason for this is well known and has been discussed in the 
Refs. J2j-^|: two possible definitions of the flow, Eckart's and Landau's definition. If we are following the 
energy flow, we satisfy exactly the energy and momentum conservation, but violate the net baryon current 
conservation. (Otherwise, if we were to choose baryon flow, we would violate the energy-momentum 
conservation.) 

The exact conservation of the energy and momentum gives for the final rapidity: 



2 (M 2 (l + eg) + 2clvl) + y/(M*(l + eg) + 2c%vlY + g§pg = „g) 
C ° &h = 2(l + eg)(M2-,g) ' (39) 



where we neglected B A If next to Q / A x A y and introduced the notation M = (1% + h)/(h — h)t 
v a = tanhj/o is the initial velocity. (The exact conservation of the baryon four-current would give: 
tanh^-B = v /M , -» cosh 2 y f , E = M 2 /(M 2 - v 2 ) ). 

It is interesting to analyze these equations, as functions of l\ and 1%. If l\ or h —> then M 2 — > 1, 
and \yf,E\ Ho- To calculate this limit for \yf t L\ we should put eg = 0, since we do not have collisions 
and, consequently, do not create QGP, thus — > yo- So there is no stopping as expected, because 
there is no reason to stop. If l\ — > h M 2 — > 00 and both expressions give Vf,E,L — * 0, i.e., complete 
stopping. So, we see that Landau's and Eckart's expressions behave similarly, have the same limits for 
minimal and maximal stopping. 

For the following part of this work we choose Landau's convention, y/ = Uf,L, which is justifiable for 
RHIC and SPS energies, where the evolution of matter is not dominated by the net baryon charge, unlike 
at lower energies where the baryon mass is still dominant and pair creation is of little importance. 
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In this case the expressions for the e / and n / are: 
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FIG. 2. The typical trajectory of the ends of two initial streaks, corresponding to numbers of nucleons, ni and 
112. Stars denote the points, where yi = yf. From t = to till those streak ends keep going according to eqs. ( J3T| , 
p8|). Later the final streak starts to move like one object with rapidity, yf, eq. (p9|) in CRF. 



The typical trajectory of the streak-ends is presented in Fig. |^. From t = t they move according to 
eqs. ( |37| , |3§| ) until they reach the rapidity y,; = yf. Later the final streak starts to move like one object 
with uniform rapidity, y/, until we reach the time when the fluid dynamical calculation starts. 

The time and position of final streak formation can be found from the condition: 



Vi = Vf 



(42) 



which gives for the ith nucleus (x\ = x , x-i = x + ) 
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(43) 
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FIG. 3. The Au+Au collision at eo ~ 65 GeV/nucl, b — 0.5(ri + rz), A — 0.08 (parameter A was introduced 
in eq. (p(j)), y = (ZX plane through the centers of nuclei), E = T 00 , laboratory frame. We note that the final 
shape of the QGP volume is a tilted disk w 45°, and the direction of the fastest expansion will deviate from both 
the beam axis and the usual transverse flow direction and might be a reason for the third flow component, as 
argued in |2a]. 



V. INITIAL CONDITIONS FOR HYDRODYNAMIC CALCULATIONS 



In this section we present the results of our calculations. We are interested in the shape of the QGP 
formed when string expansions stop and their matter is locally equilibrated. This will be the initial state 
for further hydrodynamic calculations. The time, r, at which we assume the system to reach overall local 
equilibrium and to start hydrodynamic description, is a second (after A) free parameter of our model. 
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Of course, t should be larger than the time of final streak formation, at least in the central most hot 
and dense region. For the peripheral streaks the string tension is low, and the transparency is large, 
but peripheral matter does not play a leading role in further hydrodynamic expansion. So, to have a 
homogeneous output for each streak-streak collision, we will also build the final streaks (yt, n/, e/) for 
peripheral streak-streak collisions, with lengths, A//, corresponding the lengths of the interacting region 
at the moment t = t, even if the final rapidity, yf, was not yet achieved for this particular collision. 

Au+Au, e q =65 GeV/nucI, b=0.25(r 1 +r £ ) J A=0.08 




-5 □ 5-5 5 20 40 60 80 

E, GeV/fm 3 

FIG. 4. The same as Fig. ^, but b — 0.25(ri + r^)- We see that for more central collisions the energy density 
is much larger. The QGP volume has a shape of tilted disk and may produce a third flow component [^8| . 

We may see in Figs. 0, 0, 0, that finally a QGP forms a tilted disk for 6^0. Thus, the direction 
of fastest expansion, the same as largest pressure gradient, will be in the reaction plane, but will deviate 
from both the beam axis and the usual transverse flow direction. So, the new flow component, called 
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" antiflow" or " third flow component" , will appear in addition to the usual transverse flow component 
in the reaction plane. With increasing beam energy the usual transverse flow is getting weaker, while this 
new flow component is strengthened. The mutual effect of the usual directed transverse flow and this 
new "antiflow" or "third flow component" contribute to an enhanced emission in the reaction plane. This 
was actually observed and studied earlier. One should also mention that both the standard transverse 
flow and new " antiflow" contribute to " elliptic flow" . 

Au+Au, e n =65 GeV/nucI, b=0.5(r 1 +r £ ), A=0.065 

15r t=5fm 




D 20 40 

E, GeV/fm 3 



FIG. 5. The same as Fig. [| but A — 0.065. The energy density is smaller, but the QGP volume has a similar 
shape of a tilted disk « 45° and may produce a third flow component [Eq|. We start plotting our results later 
than in Fig. H, because for smaller a the deceleration is smaller, and, so, the final streak is formed later. 



The last subplots in the Figs. [| f|, |[ || present the energy density distribution in the laboratory 
frame, E max w 50 — 90 GeV/ fm 3 for 6 = 0. It seems to be bigger than what one can expect from 
estimation based on the Bjorken model. One should, nevertheless, keep in mind that our "fireball" is not 
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homogeneous in xy plane. The average energy density for the equivalent homogeneous "fireball" would 
be lower - < E >= 22 — 29 GeV/ fm 3 . Other hydrodynamical models had to use similarly high initial 
energy density to reproduce the observed flow, e.g. in |^9| e = 23 GeV/ fm 3 has been used. 

VI. CONCLUSIONS 



Based on earlier Coherent Yang-Mills field theoretical models and introducing effective parameters, 
based on Monte-Carlo string cascade and parton cascade model results, a simplified model is introduced 
to describe the pre fluid dynamical stages of heavy ion collisions at the highest SPS energies and above. 
The model predicts limited transparency for massive heavy ions. 



Au+Au, e Q =65 GeV/nucI, b=0.25(r 1 +r 7 ) J A=0.065 
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t=3 fm 
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20 40 

E, GeV/fm 3 



FIG. 6. The same as Fig. ^, but A = 0.065. We see that for more central collisions - compared to Fig. |B|- the 
energy density is much larger, but it is smaller than in Fig. ^, because of smaller stopping. The QGP volume has 
a shape of tilted disk and may produce a third flow component |E6j|. 
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Contrary to earlier expectations — based on standard string tensions of 1 GeV/fm which lead to the 
Bjorkcn model type of initial state — effective string tensions are introduced for collisions of massive 
heavy ions. The increased string tension is a consequence of collective effects related to QGP formation. 
These collective effects in central and semi central collisions lead to an effective string tension of the order 
of 10 GeV/fm and consequently cause much less transparency than earlier estimates. The resulting initial 
locally equilibrated state of matter in semi central collisions takes a rather unusual form, which can be 
then identified by the asymmetry of the caused collective flow. Our prediction is that this special initial 
state may be the cause of the recently identified "antiflow" or "third flow component". 

Detailed fluid dynamical calculations as well as flow experiments at semi central impact parameters for 
massive heavy ions are needed at SPS and RHIC energies to connect the predicted special initial state 
with observables. 
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Au+Au, e n =65 GeV/nud, b=0.25(r+rA A=0.08, t=5 fm 
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FIG. 7. The Au+Au collision at eo = 65 GeV/nucl, t = 5 fm, A — 0.08 (parameter A was introduced in (|26|)), 
b — 0.25(ri + r?) (in our case n = r<z = r). We see that the more central plane we look at, the more nucleons 
take part in the streak-streak collisions, and therefore the more energetic and compact QGP is formed. 



Appendix A. INITIAL CONDITIONS AFTER STRING CREATION 



Our conserved quantities are 



Qo - / T 00 dV = Ax A y I T u »dz , 



=.00, 



(A.l) 



Q 3 = / T U W = Ax A y / T dz , 



(A.2) 

T 00 and T 03 are given by eq. (|3^). Before string creation the initial values of the modified energy- 
momentum tensor, T^", are 



f 1 00 -r 2 00 = eo cosh 2 yo=(-) n Q m, 

\mJ 



(A.3) 
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f 2 03 = -f? 3 = e tanh y cosh 2 y = ( — ) 2 n mv (A A) 

where m is the nucleon mass, e$ is the initial energy per nucleon, and we have used cosh 2 j/o — 7o — (m) 2 - 
vq = tanh 2/0 is the initial velocity, v = 1 is a good approximation for ultra-relativistic heavy ion collisions. 
So, 

Qo = AxAy( £ -^) 2 n m(h+l 2 ) , (A.5) 

Qz = Ax A y (—) n m(l 2 - h)v , (A. 6) 

where l\ and l 2 are the initial lengths of streaks (see Fig Q), Ax, Ay are the grid sizes in x and y directions. 
After string creation 



T = e\ cosh yi + c^ei sinh y\ + e 2 cosh y 2 + c^e 2 sinh y 2 
+ ^a 2 + B + ^n oe y° (1 - e^) + ^n e M (l - e^ 2 ) , (A.7) 



T = ei(l + Cq) cosh yi sinhj/i + ea(l + c ) cosh 2/2 sinh 7/2 



era; „, , n„, , ax 



n e yo (1 + e 2yi ) + —n e V0 (1 + e' 2y2 ) , (A.. 
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Au+Au, e n =65 GeV/nucI, b=0.5(r+rA A=0.08, t=5 fm 
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FIG. 8. The same as Fig. @, but b = 0.5(n + T2)- The stopping is smaller, consequently the QGP volume is 
less dense and less compact. 

At the point of complete penetration of streaks, t = to = (h + h)/2 (see Fig P, we introduced 
energy densities e\(to) and e2{to). We assumed transparency, i.e., that complete penetration happened 
so fast, that held, itself created during this time, did not have time to stop partons. So, the rapidities 
2/1(2) (to) = — 2/0(2/0) correspondingly, and the proper baryon densities did not change, and, thus, the 
baryon current conserved automatically. This assumption differs from how it was done in Refs. 
but they seems to be more physical, and do not change hnal results very much. Terms proportional to 
e 2yi = e~ 2y2 — e~ 2y ° <C 1 can be neglected. Then the energy and momentum conservation laws can be 
written in the form: 



Qo 

Ax Ay 



[(1 + eg) cosh 2 j/o - cl] ( ei (t )h + e 2 (t )h) + ( y + b\ (h + h) + [l\ + l\) , (A.9) 
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Ax Ay 



[(1 + eg) cosh 2 y ] (-ei(to)ii + e 2 (to)Z 2 ) - (i? - I 



(A.10) 



We neglect C 2 , close to (1 + Cq)cos1i yg i n ec L- A.£, then eqs. (A. 9, A.1C) may be solved 



ei(*o) 



UqUI 



1 



crope 



ii 



(A.ll) 



nom %- + B h+h 
T +^'(^) 2 (l+^) 2 ' 2 ^(l+cg) 



e 2 (*o) = = - * " ^ - _ b . (A.12) 




FIG. 9. The rapidity profile of streaks in the reaction plane for Au+Au collision at eo = 65 GeV/nucl, 
A — 0.065, y = 0. The rapidities of the final streaks in CRF are calculated according to eq. (|39[). Our 
profiles are in agreement with schematic sketch in paper 
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Appendix B. THE ANALYTICAL SOLUTION OF THE MODEL 

For x ± > xq we should solve eqs. (23), based on boundary conditions (p3|). Eqs. ( p3| ) leads to the 
system of equations: 

d- (h 1+e - 2yi ) + ad+h 2+ = -2an e yo + 2an e Vo e- 2y2 , (B.l) 

ad-h 1+ + d+ (h 2+ e 2y2 ) = 2an Q e yo e 2yi - 2an e y ° , (B.2) 

where a = (1 — Cq)/(1 + Cq). It is clear, that in both equations there are two terms depending on 
independent variables, so the solution will contain two undefined constants. The next step is to take eqs. 



(B.l, |B.2| ) at the values x + = Xq and x = Xq. 

h 1+ = e 2yi (ei(t )(l + c 2 )e 2ya - a 2 ( X - - x )) , (B.3) 

ad+h 2+ =c 2 + 2an e yo (e^ 2 - e - 2yo ) , (B.4) 

h 2 + = e- 2y2 (e 2 (to)(l + c 2 )e 2yo - ai(x + - x )) , (B.5) 

ad-hi+ = Cl + 2an a e yo (e 2yi - e"^ ) , (B.6) 
where we introduced new notations 

a% = ex + 4cm sinhyo , a 2 = c 2 + 4an sinh y (B.7) 

and two new constants 

ci = a(h 1+ )'\ Xa , c 2 = a(h 2 +y\x , ( B -8) 
which will be estimated by assuming a linear development for the enthalpy densities, hi+ and h 2+ , from 

t = (x± = \z(0)\) tO t = t (x± = X ). 
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c t = a{{l + c 2 Q )e t (t ) - e )/2t Q . (B.9) 
The complete analytical solution found to be 

h i+ = e(-y +1 ^ ei (t )(l + c 2 )e^ (\ ^-^j , (B.ll) 



m = n e va e { - y+1 ^ , (B.12) 
where x 1 = x~ , x 2 = x + , i, j = 1, 2 , i / j, 

6j = eta, + 2an e Vo , (B.13) 

d l ^c l - 2an e- yo , (B.14) 
ei(to)(l + cg) 



n 



e- 2 y«aj 



(B.15) 



Appendix Q. TRANSFORMATION TO THE LIGHTCONE COORDINATES 

In this Appendix we present the transformation matrixes between (t, z) and (+, — ) coordinates. Indexes 

k run for 0, 3, indexes a, (3, 7 are used for +, — . 
The transformation of the coordinate system is 

x ± =t±z, (C.l) 

thus, for all contravariant vectors we have the same: 

V ± = V° ± V 3 . (C.2) 

For (t, z) coordinates we have 
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Then 



so, 



g ik =g lk =\ _ x ■ (C3) 



<m = ^=(11), (c.4) 



J _ „.a _ 1 f 1 1 



2 V 1 1 



2 

The backward transformation is 



so, 



(C.5) 



9 a i=9i=[ \ _\ J ■ (C6) 



T a/3 = g^gf , (C.7) 



T++ = 1 (T 00 + T 03 + T 30 + T 33 ) , (C.8) 



yH — _ ^_ (rpQO _ rp03 _|_ y30 _ ji33~J ^ g-j 



T -+ = 1 (yOO + ^03 _ T 30 _ T 33^ ^ ( C 1Q ) 



ji— _ 1 ^y 00 — T 03 — T 30 + T 33 ^j (C 11) 



T ij = gi a T°<P g e, (C.12) 



T 00 = i (T++ + T+- +T-+ +T— ) , (C.13) 



T 03 = 1 ^ r ++ _ T +- + T -+ _ T — ■) ^ ( C U ) 
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T 



30 



T 4 



. — 'J 1 — r r £~ 



J~<33 _ + + r p-\~— r-p — -\- _|_ fj^- 



(C.15) 
(C.16) 
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